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V' 

The  discrete  least  squares  method  is  convenient  for  computing  polynomial 
approximations  to  functions.  We  investigate  the  possibility  of  using  this 
method  to  obtain  polynomial  approximants  good  in  the  uniform  norm,  and 
describe  applications  both  to  the  case  when  the  function  to  be  approximated  is 
known  on  a  discrete  point  set  only  and  to  the  case  when  we  can  freely  choose 
the  set  of  least  squares  nodes.  Numerical  examples  are  presented. 
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SIGNIFICANCE  AND  EXPLANATION 


The  simplicity  of  computing  least  squares  polynomial  approximations  to 
functions,  using  function  values  on  a  discrete  point  set  only,  makes  it 
attractive  to  use  the  least  squares  method  for  computing  polynomial 
approximations,  which  are  good  also  when  the  error  is  measured  in  the  uniform 
norm.  He  describe  how  this  can  be  done,  and  consider  the  cases 

1)  the  function  to  be  approximated  is  known  on  a  discrete  point  set 
only.  For  example,  we  may  wish  to  approximate  a  function  on 
[-1,1],  but  the  function  is  only  known  on  an  equidistant  point  set 
on  the  interval. 

2)  He  are  free  to  select  an  arbitrary  set  of  discrete  least  squares 
nodes.  For  example,  we  may  wish  to  approximate  an  analytic  function 
on  a  bounded  region  in  the  complex  plane,  and  the  function  is  known 
on  the  boundary  of  the  region.  He  discuss  how  to  allocate  least 
squares  nodes  on  the  boundary  without  computing  conformal  mappings. 
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ON  POLYNOMIAL  APPROXIMATION  IN  THE  UNIFORM  NORM 
BY  THE  DISCRETE  LEAST  SQUARES  METHOD 

Lothar  Reichel 

t.  Introduction 

Polynomial  approximations  to  functions  are  conveniently  computed  by  the 
discrete  least  squares  method.  He  describe  how  this  method  can  be  used  to 
compute  polynomials,  which  provide  good  approximations  also  when  the 
approximation  error  is  measured  in  the  uniform  norm  on  the  whole  domain  of  the 
function.  Important  is  the  choice  of  an  appropriate  relation  between  the 
number  of  discrete  function  values  to  be  used,  and  the  degree  of  the 
polynomial.  This  choice  should  depend  on  the  distribution  of  the  least 
squares  nodes. 

Ex.  1.1.  Consider  the  approximation  of  f(x)  s*  (1  +  25x  )  on  [-1,1]. 

k-1 

Use  function  values  at  m  equidistant  points  x^  :*  -1  +  2  ^7^,  k 
■  1(1 )m,  to  compute  a  least  squares  polynomial  of  degree  <  n. 

a)  n  -  m.  The  uniform  norm  approximation  error  increases  with 
m,  see  Runge  [5],  or  Dahlquist-Bjorck  [4],  section  4.3.4. 

b)  n  m  —  /m.  The  computed  sequence  of  polynomials  for 

✓2 

increasing  n  converges  maximally  to  f(x)  in  the  sense  of 
Halsh  [7],  i.e.  the  exponential  rate  of  convergence  is 
optimal.  Details  are  given  in  section  2. 
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In  section  2,  we  state  results  on  the  rate  of  convergence  in  the  uniform 
norm  of  the  polynomial  approximant  computed  by  the  least  squares  method.  Our 
approach  is  to  first  give  error  bounds  containing  certain  interpolation 
operator  norms,  and  then  bound  the  interpolation  operators.  Numerical 
examples  are  presented  in  section  3. 


2.  Estimates  for  the  least 


rator 


Let  T  be  a  Jordan  curve  or  Jordan  arc  in  the  complex  plane.  Denote  the 
open  set  bounded  by  T  by  Q.  If  F  is  a  Jordan  arc  then  ft  is  void.  We 
wish  to  coelute  polynomials,  which  provide  good  approximations  to  functions 


f  on  T  in  the  uniform  norm, 


(2.1) 


Ifl  sup  |f(s)| 

u  *er 


Let  {z,  }*  be  a  set  of  distinct  nodes  on  I*.  Introduce  the  inner-product 

X* I 


and  corresponding  semi-norm 


(2.2) 


(2.3) 


(f,9)  £  I  9<*kfm>  . 


Ifl  t-  (f,f) 


where  the  bar  denotes  complex  conjugation.  For  a  given  function  f  on  T, 
let  Ln^mf  denote  the  best  polynomial  approximant  of  degree  <  n  with 
respect  to  the  semi-norm  (2.3),  m  >  n.  Let  Inf  be  the  polynomial  of 
degree  <  n  interpolating  f  at  n  distinct  points  {w^}”^  on  I*.  We 


introduce  the  notation 


(2.4) 


I  -4  L 
n  n 


,m  if  {wk,k.1C{VA-1'W)t  distinct  . 


We  also  need  the  definition 


(2.5) 


E  (f)  :■  sup  If  -  p  I  , 
n  n  u 

Pn 


where  the  supremum  is  taken  over  all  polynomials  of  degree  <  n. 

Theorem  2. 1 

Define  L  and  I_  on  the  set  of  functions  continuous  on  I*  U  ft  and 

it*  n  n 

analytic  on  0,  and  let  both  domain  and  range  of  Ln  m  and  In  be  equipped 
with  the  uniform  norm  (2.1).  (If  r  is  an  arc,  then  0  is  empty. )  Then 


(2.6) 


IL  I  <  II  1(1  +  /m)  V  I  -*  L  ,  V  m  >  n  . 
n,m  n  n  n,ra 


The  growth  rate  of  the  right  hand  side  is  the  smallest  possible. 

Proof.  Let  PRf  be  a  polynomial  of  degree  <  n  such  that  IPnf  -  fl( 


E„(f).  Then 
n 

(2.7) 


IL  f  -  fl  <  IP  f  -  fl  UP  f  -  fl  “  E  (f)  . 

n,m  m  n  m  n  u  n 


Let  ML.  -  and  let  t.{z)  denote  the  Lagrange  polynomials  associated 
n  n,m  X 


with  In,  i.e. 


(2.8) 


II 

(I  f )  (z )  -  X  f(w.  )JL  (z) 
n  km1  x  K 


Express  Ln  m  with  the  same  polynomial  basis , 


(2.9) 


n 

(L  f)(z)  -  X  <VVZ> 
nrm  .  -  k  k 


for  some  constants  a^.  Substituting  (2.9)  into  (2.7)  yields 
lo^  -  f(wk)|  <  /m  En(f)  ,  )c  -  1  ( 1  )n  . 


Substitution  into 


|  (L  f )  ( z )  |  <  X  la  -  f(w.)|  |i.(z)|  +  X  lf(wk>l  iy*>l 
n,m  k-1  *  k  *  k-1  K  K 

yields 

(2.10)  IL  I  -  sup  IL  fl  <  II  I  +  II  l/m  sup  E  (f)  . 

n,m  Ifl  *1  n,m  n  n  n  |f|  n 

u  u 


This  proves  (2.6).  We  next  give  an  example,  where  the  growth  rate  of  the 
bound  (2.6)  for  increasing  m,  and  n  held  fixed,  is  obtained. 

Let  fg(x),  x  6  [-1,1],  q  >  0,  denote  the  piece-wise  linear  function 


(2.11) 


f  (x) 

q 


q>t  ,  |x|  <  Vq 

x/lxl  .  Vq  <  IjcI  <  1 


Approximate  fg  by  a  1st  degree  polynomial  L2,mfq  on  [-1,1]  for  m 
even.  The  polynomial  L-  f  we  define  by  least  squares  approximation  at 

“I*® 

1  1  n 

the  m  points  x1  «  ~,  x2  -  x3  -...■  x^  •  -,  x^^  “  ~xy  j  *  HD 

/2  2 


Then 


1  +  f-  -  1 ) 1*- 

(L-  f  )(x)  - - - - r-3  x 

2'“  q  1  +2(t-  1)»^2 

2  q 


Letting  q  “  /  —  -  1 ,  we  obtain 


(L-  f  )(x)  -  -(1  +  /m-1 ) x  . 

<ii>  q  •> 

For  m  sufficiently  large, 

(2.12)  IL  I  >  IL  f  I  >  IL  f  -  f  I  -Ifl  *  -^Vjn-T  -  2)  -  1  . 

n,m  n,m  q  u  n,m  q  q  u  q  u  3 

Moreover, 

x-x1  x-x  , 

(2.13)  inf  1X_I  -  l| - —|  +  1 - —|  I  <  z  1°(2)  +4  V  m  >  4  . 

.  ,,  2  'x  -x,  *x4-x  u  v 

I3L,  ml  1m 

2  2  fin 

Combining  (2.12)  and  (2.13),  we  see  there  is  an  I2  and  a  constant  d  >  0, 


such  that 


(2.14) 


IL.  I  >  II  l/m  d  V  a  sufficiently  large  . 
2  ,m  2 


The  least  squares  nodes  of  this  example  are  not  distinct,  but  we  can  find 

distinct  nodes  close  to  those  used  in  this  example  and  such  that 

IL-  I  >  IX  I/m  •  d/2. 

2fin  2 


Generally,  we  will  select  n  as  an  increasing  function  of  m. 
Additional  smoothness  of  the  functions  to  be  approximated  will  decrease  the 
growth  of  IL  I  with  m,  n  ■  n(m).  The  next  theorem  illustrates  this. 

flfB 

Theorem  2.2 


a)  Let  T  “  t— 1 ,1J  and  let  L_  _  have  the  domain  F.  ,  :■ 

n,m  d,k 

{f  j  f  e  Ck[-1,1],  1^— I  <  d}  equipped  with  norm  (2.1).  Then,  for  some 

dz  u 

constant  D  depending  on  the  constant  d  and  integer  k, 

IL  I  <  IX  1(1  +  0m1/2n"k)  V  I  4  L 
n,m  n  n  n,ro 


(2.15) 


b)  Let  r  be  a  smooth  Jordan  curve  with  interior  Q.  Define  L_  _  on 

n^m 

V 

Fj  l  *  (f  i  f  analytic  in  G,  f  6  c  (G  U  D,  I — • 1  I  <  d}  equipped  with 

_  K  U 

dz 

norm  (2.1).  Then,  there  is  a  constant  D  such  that 

(2.16)  IL  I  <  II  1(1  +  D/m  ,  v  I  -(L  ,  n  >  k+2  . 

n,m  n  n  n  n,m 

Proof .  The  inequalities  (2.15)  and  (2.16)  follow  from  (2.10).  To  obtain 
(2.15)  we  use  Jackson's  theorem,  see  Cheney  [2],  p.  147.  (2.16)  follows  from 

Smirnov-Lebedev  [6] ,  p.  99. 

In  the  following  we  will  assume  that  there  is  a  distribution  function 
S(z)  on  T  with  a  strictly  positive  derivative  for  the  least  squares  nodes 
z ^  ,  and  that  these  nodes  are  defined  by 

(2.17)  S(z.  )  =  k  ■  1(1)m,  m  ■  1,2,... 

K,m  m 

The  sets  of  interpolation  points  will  be  subsets  of  the  least  squares  node 
sets. 

The  next  theorem  gives  estimates  for  II  I  for  various  distributions  of 

n 

interpolation  points.  We  first  single  out  two  special  point  distributions. 
The  definitions  follow  Walsh  [8] ,  ch.  7. 

Definitions 

Let  G(z)  be  the  Green's  function  for  the  region  exterior  to  V  with  a 

3 

logarithmic  singularity  at  infinity.  Let  t~*  denote  the  outward  normal 

on 

derivative.  A  point  set  on  T  is  called  uniformly  distributed  if 

the  8  defined  by 
X  $  n 

,2-,a)  £  “’i*'1  -  v.  •  k  • 

for  any  0  <  c  <  d  <  1,  satisfy 

lim  —  {Number  of  9. 
n  k,n 


in  [c,dj }  ■  d-c 


Integration  in  (2.17)  is  understood  to  be  carried  out  along  r  in  the 
positive  sense. 

It  8  -  Vk,  on  8  -  K  Vk,  then  the  point  set  {w  }"  .  is 

K|R  n  K#n  n  k  k* i 

said  to  be  equidistrlbuted  on  I*. 

Ex.  2.1.  The  zeros  of  the  Chebyshev  polynomial  TR(x)  -  cos(n  arc  cos(x)) 
are  equidistrlbuted  on  [-1*1].  Equidistant  nodes  on 
r  -  (s  : | z|  ■  1}  are  equidistrlbuted. 

Theorem  2.3 

If  the  interpolation  points  are  equidistrlbuted,  and  I*  is  an  interval, 
or  an  analytic  Jordan  curve,  then  there  is  a  constant  3  such  that 

(2.19)  II  I  <7  ln(n)  +6,  n  -  0,1,2,...  . 

n  w 

If  the  interpolation  points  are  uniformly  distributed  and  T  is  a  Jordan 
curve  or  Jordan  arc,  then 

(2.20)  II  l1/n  ♦  1,  u  ♦  •  . 

n 

If  the  interpolation  points  w^  »  wkn,  k  *  1(1 )n,  in  the  limit  n  ♦  •  are 
distributed  according  to  a  density  function  o,  then 

(2.21)  II  |1/n  -  ~  ,  n  ♦  -  , 

"  P0 

where  P1  :«  exp (sup  /  ln|z-?| o( O | d?j ) ,  and  pQ  :* 

*er 

expdnf  L  ln|s-(|o(()|dd). 

*er 

Proof. 

Statement  (2.19)  is  well-known  for  T  being  an  interval.  For  T  being 
an  analytic  Jordan  curve,  (2.19)  follows  from  Curtis  [3],  theorem  B.  (2.20) 
is  a  special  case  of  (2.21).  The  latter  can  be  shown  by  potential  theoretic 
methods.  This  is  carried  out  in  the  appendix. 


Approximation  method 

Guided  by  theorem  2.3,  we  select  n  -  n(m),  so  that  among  the  least 

squares  nodes  {z.  }™  there  is  a  subset  of  n  points,  which  is  uniformly 

kfn  k* i 

distributed  on  T  as  m  •.  The  assumed  existence  of  a  distribution 
function  S(z),  with  strictly  positive  derivative,  for  the  z^  m  guarantees, 
that  such  a*  function  n(m)  does  exist,  and,  moreover,  we  can  let  n(m)  ♦  « 
as  m  ♦  •. 

Maximal  convergence,  defined  in  ex.  1.1,  is  obtained  if  IL  l^n  ♦  1 

nfn 

as  n,  m  ♦  •».  From  (2.6),  we  see  that  a  sufficient  condition  for  maximal 
convergence  is  that  the  function  n(m)  is  such  that  n(m)  >  ma  for  some 
constant  a  >  0.  The  next  theorem  provides  an  example. 

□ 

Theorem  2.4 

Let  T  *  [-1,1],  and  let  z.  ■  -1  +  2~  ",  k  m  1{1)m.  Then  n  m  —  /m 

x,m  m 

satisfies  the  proposition.  Conversely,  there  is  no  subset  of  n  **  cm 
uniformly  distributed  nodes  as  m  ♦  <*,  for  any  constant  c  >  0. 

Proof.  See  appendix. 

□ 

Remark 

An  approximation  problem  related  to  that  treated  in  Theorem  2.4  has  been 
discussed  by  Bjork  [1],  who  measured  the  approximation  error  I^^f  -  f  with 
a  semi-norm  1*1^,  with  £>>  m.  The  computation  of  the  Euclidean  norm  of 
matrices  defining  the  change  between  different  orthonormal  polynomial  bases 
on  [-1,1],  led  Bjork  to  the  suggestion  that  n  should  be  selected  4  2/m. 
This  choice  of  n  is  close  to  ours,  since  —  «  2.22. 


□ 


1^ 


3.  Numerical  examples 

When  applying  the  approximation  method,  we  discern  the  following 
different  cases. 

a)  the  function  f  is  known  on  all  of  T. 

b)  the  function  f  is  known  only  on  a  finite  point  set  on  I*. 

c)  the  normal  derivative  of  the  exterior  Green's  function  for  T  is 
explicitly  known 

d)  the  normal  derivative  of  the  exterior  Green's  function  for  r  is  not 
simply  available. 

Only  if  a)  and  c)  are  true,  or,  if  we  have  the  situation  described  in  theorem 

2.4,  then  we  know  how  to  select  n  “  n(m).  In  general  we  compute  LR  mf  for 

several  n  <  m.  This  is  illustrated  in  example  3.3.  If  then  a)  holds  true, 

we  can  compute  IL  f  -  fl  and  select  the  best  of  the  computed  L  _f’s. 

n,m  u  n,m 

If  instead  b)  is  true,  the  selection  of  an  appropriate  Ln  mf  must  be  based 
on  a  numerical  perturbation  analysis. 

Our  first  example  is  a  continuation  of  examples  1.1  and  1.2.  All 
computations  were  carried  out  on  a  AX/780  in  double  precision,  i.e.  with  15 
significant  digits. 

Ex.  3.1.  Let  r  -  [-1,1],  f(z)  -  (1  +  25z2r1,  and  z.  m  -  -1  + 

f  •p 

k  +  1(1 )m.  Let  n  be  the  largest  integer  <  —  /m.  I  f  below 

z  n 

denotes  the  interpolation  polynom  Ifif  determined  by 

interpolating  f  at  the  zeros  of  Tn(x)  *  cos(n  arc  cos  (x)). 

T 

I^f  is  known  to  converge  maximally  to  f.  In  all  examples  on 
approximation  on  [-1,1],  we  have  used  the  polynomials  T^(x)  as 


-9- 


basis  functions 


n 


m 


il  e  -  f» 
n,  m  u 


IL 


f  -  fl 


1/n 


uTf  -  nv° 

n  u 


10 

7 

0.688 

1.15 

40 

14 

0.173 

1.02 

160 

28 

0.161 

•  10-1 

1.03 

640 

56 

0.904 

•  io“4 

1.02 

2560 

112 

0.192 

•  10"8 

1.01 

The  entries  of  the  last  column  are  close  to  1  and  decrease  as  n,  m 
increase.  This  indicates  the  maximal  convergence  of  Ln<mf  to  f. 


□ 


Ex.  3.2.  We  continue  the  previous  example.  Increase  the  value  of  n  to 
n  *  0.5m,  everything  else  kept  as  in  ex.  3.1. 


m 

n 

IL  f  - 

n,m 

10 

5 

0.319 

40 

20 

0.523 

160 

80 

12.221 

The  error  increases  exponentially. 

Ex.  3.3.  Let  fl  be  the  interior  of  the  curve  T  s*  {x+iy,  x(t)  :* 
y  cos(t)  +  y  cos(2t)  -  y,  y(t)  sin(t)  +  j  sin(2t)  - 
~  sin(4t) ,  0  <  t  <  2ir} . 


□ 


-10- 


We  wish  to  approximate  functions 


analytic  on  0  and  continuous  on 
(1  U  T  by  polynomials.  T  is  some¬ 
what  similar  to  the  ellipse  &  :« 

{*,  s(t)  *  x(0)cos(t)  +  iyt'jjsin  (t) 
0  <  t  <  2*}. 


The  point  set  {stt^  i*  «<iuidi®tributed  on  &  if  t^  B  “ 

k-1 

2ir  k  “  1  ( 1 )m.  This  leads  us  to  select  the  least  squares  nodes  m  “ 

k-1  k-1 

x(2w  - )  +  iy(2w  - ),  k  *  1(1  )m,  on  T.  As  basis  functions  we  use  the 

IB  IB 

Chebyshev  polynomials  of  the  1st  kind  for  the  interval  between  the  foci 
of  &  and  we  scale  them  so  that  their  maximum  magnitude  on  T  is  1.  This 
basis  is  sufficiently  well-conditioned  to  allow  representation  of  polynomials 
of  a  fairly  high  degree.  Let  f(z)  *  (z  -  2x{*))  .  For  m  *  80,  we 

compute  Ln^m  for  n  **  30(10)70.  The  best  approximation  so  obtained  is 
underlined. 


m 

n 

IL  f  -  fl 
n,m  u 

80 

30 

0.193  •  10"3 

80 

40 

0.914  •  10“5 

80 

50 

0.437  •  10'6 

80 

60 

0.250  •  10”5 

80 

70 

0.241  •  10"2 

11- 


How  overdetermined  the  linear  system  should  be  to  yield  the  smallest 


approximation  error  depends  on  the  location  of  the  singularities  of  the 
function  to  be  approximated.  Reflecting  the  singularity  of  f  in  the 
imaginary  axis,  we  obtain  g(s)  ■  (s  -  2x(0))-1. 


m 

n 

IL  g  -  gl 

n,m’  3  u 

40 

20 

0.413  •  10"3 

40 

30 

0.114  •  10"4 

40 

39 

0.611  •  10~® 

40 

40 

0.592  •  10”6  . 

□ 
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Appendix 


Proof  of  theorem  2.3. 

By  potential  theoretic  methods,  we  prove  (2.21)  of  theorem  2.3.  Let  the 
interpolation  nodes  be  and  consider  the  Lebesgue  function 

Q 

W. 


(A.  1 ) 


A  (z)  «' 
n 


n  n 

l  « 

k-1  j-1  ,wk  wj 

i+k 


For  z  +  Wj,  we  have 


(A.2 ) 

1  n 

1  In  n 

n  j-1 

Ir-Wjl 

i  n 

"I  l 

j-1 

In | z-w .  |  = 
3 

/rm|z  - 

w |  o(w)  |  dw|  +  0(ln^}-)  > 

j+* 

and 

(A. 3 ) 

-  In 
n 

n 

B  1 

j-1 

V^j1  * 

/r  ln|wk-w 

|  cr(v)  |dw| 

+0(lniHl)  . 

n 

j+* 

Substituting  (A. 2),  (A. 3)  into  (A.1)  yields 

A  (z)  «  n2  exp(n  J_  ln|z-w|o<w) |dw| )  •  0(1)  • 
n  * 


j  exp(-n  Jj,  ln|wk-v|o(w)|dw|),  z  e  T  . 


Hence, 


and  finally. 


2  n 


-n 


A  (z)  ■  n  P1  •  0(1)  •  n  P0  , 


II  lV"  -  IAl’/n  ■  ♦ 

n  u  pQ 


If  the  interpolation  nodes  are  uniformly  distributed  as  n  ♦  •»  then  c(w) 
is  the  equilibrium  density  function  for  T,  and  in  this  case  “  PQ. 


Proof  of  theorem  2.4. 

The  zeros  of  the  Chebyshev  polynomial  Tn(x)  -  cos(n  arc  cos(x))  are 
equidistributed,  and  therefore  uniformly  distributed,  on  [-1,1].  Let 
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NT(c,d)  denote  the  number  of  zeros  of  T  (x)  in  the  interval  tc,d] , 
n  n 

-1<c<d<1.  Given  a  set  of  n  uniformly  distributed  nodes  on  [-1,1], 

let  N„(c,d)  be  their  number  in  [c,d],  -1  <  c  <  d  <  1.  Then 
n 

lim  -  N  (c,d)  -  lim  -  NT(c,d). 

n  n  ^  n  n 

n+«»  n+« 

We  first  show  that  n  cannot  be  equal  to  cm,  for  any  constant  c  >  0. 

Assume  the  contrary,  i.e.  there  is  a  constant  c  >  0,  such  that  as  n  +  «•, 

there  is  a  subset  of  n  *  cm  uniformly  distributed  nodes.  For  any  e  with 

0  <  e  <  1,  we  then  have 

lim  -  (N  (-1,-1+e)  +  N  1 1-e,  1 ) )  <  e  , 
in  n  n 

and 


lim 

n+« 


n 


N  (-1+e,1-e) 


2 

z 


arc  sin(l-e) 


Therefore 

lim  (— (N  (-1+e, 1-e)  +  N  (-1,-1+e)  ♦  n  ( 1-e, 1 ) )  <  -  +  -  arc  sin(l-e)  . 

„+«  n  n  n  n  c  * 

It  remains  to  be  shown  that  to  each  c  >  0,  there  is  an  e,  0  <  e  <  1,  such 

that  —  +  z  arc  sin(l-e)  <  1.  Let  h.(e)  :•  and  h_(e)  :* 

C  V  1  C  A 

1  -  —  arc  sin(l-e).  Then  h.(0)  *  h,(0)  »  0.  Further  h'(e)  “  — , 

IT  1  Z  1C 

2  1 

h*  (e)  ■ - ,  e  >  0.  In  a  punctured  neighborhood  of  e  •  0,  we  have 

2  i  -  - 

/2e-e 

h^e)  >  h^(e).  This  shows  that  h2(e)  >  h^(e)  in  that  neighborhood,  and  we 
are  ready. 

Next  we  show  that  a  subset  of  —  /m  nodes  can  be  distributed  uniformly 

a 

dG  1  1 

with  respect  to  t-  (z)  - - - ,  z  e  [-1,1].  The  distance  between  the  last 

,n 


two  largest  zeros  x1>n,  x2,n  °*  Tn^x)  i® 


2 

x1,n  ‘  X2,n  “  co,(fe)  -  COB(S)  *  2  sin(J,Bin(fe)  "  h  +  0<~4)  * 


With  m  equidistant  nodes#  z  •  -1  ♦  - #  k  *  1(1 )m#  we  have 

k#m  n 


(A.  4) 


2  * 


2  ' 


or  equivalently  n  »  —  a  .  Since  the  zeros  of  T_(x)  are  most  dense  at  the 

VS 

ends  of  the  interval,  the  choice  (h.4)  of  n  guarantees  that  for  every  m  a 

subset  of  n  nodes  can  be  selected  which  is  close  to  the  set  of  seros  of 

dG 

Tn(x),  and  uniformly  distributed  with  respect  to  ■gj  as  n  ♦  •. 
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